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摘要 ， 针对 目前 电压 暂 降 检测 时 存在 的 由 于 噪声 和 谐 波 等 干扰 导致 的 分 析 精 度 差 、 
实时 性 不 高 等 问题 ， 提 出 一 种 基于 双 小 波 变 换 的 电压 暂 降 检测 方法 。 该 方法 首先 采用 
db20 小 波 对 信号 进行 多 尺度 分 解 ， 提 取 含 暂 降 信息 的 基 波 信号 。 然 后 采用 db10 小 波 
对 基 波 信号 进行 单 尺度 分 解 ， 检 测 高 频 分 量 的 模 极 大 值 对 电压 暂 降 实现 定位 ， 从 而 降 
低 了 噪声 和 谐 波 等 干扰 信号 的 影响 ， 提 高 了 小 波 分 析 的 检测 性 能 。 
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Study of Voltage Sag Detection Method Based on 
Double Wavelet with Anti-Interference Feature 


XuDan Huang Lijun Yang Xiaohui LiXianwei Wang Yi LiBo 
( Xuji Group Corporation Xuchang 461000 China ) 


Abstract: In view of the poor analytical accuracy and weak real-time performance 
caused by signal noise and harmonic wave exist in the voltage sag detection, a 
comprehensive voltage sag detection method is proposed to weaken the impact of the 
interference signals by the double wavelet analysis. This method decomposes the signal 
at multi-scale with db20 wavelet to catch the fundamental wave, and then decomposes the 
fundamental wave at single-scale with db10 wavelet to detect the modulus maximum of 
high-frequency components. The proposed method can suppress the interference signal 
and improve the performance of voltage sag detection by wavelet analysis. 
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1 引言 


近年 来 ， 随 着 电力 电子 技术 的 快速 发 展 以 及 
计算 机 和 自动 化 控制 设备 在 工业 和 社会 领域 的 大 规 
模 应 用 ， 人 们 对 供电 系统 的 电能 质量 提出 了 更 高 要 
求 。 电 压 暂 降 是 电能 质量 中 的 主要 问题 。 根 据 我 国 
国家 电能 质量 标准 ， 电 压 暂 降 是 指 在 工 频 条 件 下 ， 
电压 方 均 根 值 突然 降 至 0.1 ~ 0.9pu， 并 在 短暂 持续 
10ms ~ lmin 后 恢复 正常 的 现象 中。 电压 暂 降 对 敏 
感 负荷 造成 了 极 大 危害 ， 如 引起 PLC 失灵 ， 计 算 机 
存储 数据 丢失 等 ， 导 致 巨大 经 济 损失 " ， 因 而 受到 
广泛 关注 。 

目前 ， 检 测 电压 暂 降 信号 的 方法 主要 有 时 域 法 
和 坐标 变换 法 。 时 域 法 通过 求 取 差 变 信号 检测 电压 
暂 降 发 生 的 时 刻 及 持续 时 间 ， 坐标 变换 法 包括 短 时 
傅 里 叶 变 换 、dg 坐标 变换 以 及 小 波 变 换 等 中"。 其 
中 ， 小 波 变换 具有 恨 好 的 暂 态 信号 分 析 能 力 ， 能 
够 准确 识别 电压 突变 点 的 特征 量 ， 受 到 研究 人 员 
的 广泛 关注 。 文 献 [12] 采用 现场 可 编程 序 阵列 
(Field-Programmable Gate Array，FPGA) 实现 了 
小 波 变 换算 法 对 电压 暂 降 的 实时 检测 ， 时 间 延 迟 达 
到 lms， 满 足 工程 应 用 的 实时 性 要 求 。 然 而 ， 滤 波 
和 去 噪 是 确保 电压 暂 降 检测 实时 和 精确 的 关键 。 在 
实际 应 用 中 ， 检 测 信 号 中 存在 的 噪声 和 谐 波 干扰 会 
影响 小 波 变 换 的 检测 分 析 能 力 ， 导 致 误差 增 大 其 至 
造成 误 判 。 文 献 [13] 提出 了 采用 设置 闵 值 的 方法 对 
小 波 高 频 系数 进行 量化 以 消除 噪声 影响 ， 而 国 值 的 
选取 在 实际 应 用 中 往往 难以 确定 。 文 献 [14] 提出 了 
运用 多 尺度 小 波 分 解 的 方法 ， 充 分 考虑 利用 高 频 和 
次 高 频 分 量 来 抑制 噪声 幅 值 ， 但 对 于 较 大 噪声 幅 值 
以 及 存在 谐 波 的 情况 下 容易 发 和 后 误 判 。 文 献 [15] 提 
出 了 利用 加 权 最 小 二 乘 估计 算法 和 重 置 协 方差 来 快 
速 判 断 电压 暂 降 ， 但 谐 波 分量 的 频繁 变化 将 会 导致 
国 值 不 断 重 置 。 文 献 [16] 提出 了 利用 双 小 波 进 行 粗 
估计 和 精确 定位 起 止 时 刻 的 方法 ， 以 排除 噪声 干扰 ， 
但 高 频次 谐 波 存在 影响 模 极 大 值 的 奇异 点 ， 降 低 了 
小 波 分 析 的 检测 能 

本 文 首 先 研 究 了 噪声 和 谐 波 在 小 波 各 尺度 上 的 
传播 特性 ， 提 出 了 双 小 波 检测 思想 ， 即 在 大 尺度 下 
分 析 提 取信 号 基 波 成 分 ， 在 小 尺度 下 分 析 电 压 突变 
信息 。 在 此 基础 上 ， 本 文 研究 了 消失 和 矩 阶 数 对 基 波 
信息 提取 精度 和 和 暂 降 信息 定位 的 影响 ,选取 合理 的 
小 波 函 数 ， 确 定 双 小 波 检测 方案 ， 并 通过 Matlab/ 
Simulink 仿真 验证 了 该 方法 的 可 行 性 。 
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2 小波 分 析 方 法 


2.1 小 波 变换 

小 波 变换 的 窗口 大 小 可 根据 信号 频率 进行 调整 ， 
其 本 质 思想 是 将 信号 表示 为 以 小 波 函 数 和 派生 出 的 
尺度 函数 为 基 的 代数 和 。 对 于 任意 信号 f(D)EL(R) 
和 小 波 函 数 到 四， 其 连续 小 波 变 换 可 表示 为 


WT.(a,b)=< f(t), 00>- 工 | ， /ov [Sh (1) 


式 中 ，< f(D)，Wis(D) > 表示 f(D 与 sl) 的 内 积 ; 
到 为 小 波 函 数 炎 的 共 轿 函数 ，a 为 尺度 因子 ， 表 征 
信号 在 频 域 的 伸缩 ， 5 为 平移 因子 ， 表 征 信号 在 时 
域 的 平移 。 
2.2 多 分 辨 分 析 与 滤波 

小 波 函 数 的 时 频 窗口 中 心 和 半径 随 尺度 因子 的 
变化 而 变化 。 通 常 ， 将 尺度 因子 进行 离散 化 (a= 
2,j 为 整数 )， 实 现 频 域 的 逐 层 下 行 二 分 。 对 信号 
进行 小 波 变 换 相 当 于 将 信号 分 解 到 不 同 频 段 中 ， 小 
波 系数 表征 各 频段 信号 的 分 量 信息 。 各 频段 的 小 波 
系数 可 以 经 过 逆 变 换 进 行 重 构 以 还 原 信号 分 量 ， 进 
而 获得 信号 的 基 波 或 谐 波 信息 。 图 1 为 尺度 为 3 (j= 
3) 的 小 波 分 解 示 意图 。 信 号 5 经 小 波 变换 后 表示 为 
S=4+Di+D+D， 4、4、43 表 征 信号 的 低频 轮 
廊 ， 对 应 频段 0 ~ /4、0 ~ /8、0 ~ 7N16， 其 中 大 
为 采样 频率 ，D,、D,、D; 表征 信号 的 高 频 细 节 ， 对 
应 频段 /4 ~ f/2、fY8 ~ f/4、f116 ~ f/8， 对 于 信 
号 的 高 频 分 量 ， 采 用 小 尺度 短视 窗 进 行 分 析 ， 对 于 
信号 低频 分 量 ， 采 用 大 尺度 宽 视 窗 进 行 分 析 。 通 过 
合理 地 设置 采样 频率 与 分 解 尺 度 ， 可 使 信号 的 基 波 
与 谐 波 有 效 分 离 ， 基 波 信息 存在 于 最 低频 段 中 ， 谐 
波 信息 分 布 于 各 高 频频 段 中 。 


可 
4 D 第 1 层 | 
a +- Se 
五 | 第 2 层 
4; D, 第 3 层 | 


1 小波 变换 分 解 示意 图 
Fig.1 The schematic diagram of wavelet transform 
2.3 奇异 性 检测 与 暂 降 定位 


电压 暂 降 属于 奇异 信号 的 一 种 ， 相 当 于 在 信和 号 
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的 奇异 点 处 琶 加 脉冲 函数 或 阶 跃 函数 。Mallat 证 明 
了 信号 的 局 部 奇异 性 与 小 波 变换 模 极 大 值 密切 相关 。 
信号 的 奇异 点 与 小 波 系数 的 模 极 大 值 点 在 时 域 上 存 
在 对 应 关系 ， 并 从 一 个 尺度 传播 到 另 一 个 尺度 ， 同 
时 在 各 尺度 上 均 保 持 极 大 值 ， 通 过 检测 小 波 分 解 高 
频段 的 模 极 大 值 点 ， 进 而 确定 原始 信号 中 暂 降 发 生 
的 起 止 时 刻 。 图 2 为 电压 暂 降 信号 的 单 尺 度 分 解 结 
果 ，D) 为 高 频段 的 重 构 信号 分 量 。 从 高 频 分 量 DD 
中 可 以 明显 识别 出 模 极 大 值 点 ， 分 别 对 应 电压 暂 降 
的 起 止 时 刻 。 


1 


图 2 小 波 变换 暂 降 定位 原理 
Fig.2 The principle of voltage Sag detection by DWT 


3 ” 双 小 波 暂 降 检 测 


小 波 分 析 方 法 通过 检测 最 高 频段 D 或 次 高 频 
段 D, 的 模 极 大 值 来 实现 电压 暂 降 定位 。 在 实际 应 用 
中 ， 由 于 噪声 和 谐 波 的 存在 ， 往 往 会 导致 检测 性 能 
下 降 其 至 失效 。 
3.1 传播 特性 

暂 降 信号 、 品 声 与 谐 波 在 小 波 各 尺度 的 传播 特 
性 各 不 相同 。 以 220V 含有 和 暂 降 的 工 频 信号 为 例 ， 
加 入 3 次 谐 波 、13 次 谐 波 以 及 白 噪 声 干扰 ， 分 别 研 
究 暂 降 、 噪 声 与 谐 波 在 小 波 各 尺度 的 传播 特性 ， 信 
号 特征 及 分 析 条 件 见 下 表 。 

表 ”信号 特征 参数 及 仿真 分 析 条 件 


Tab. Signal characteristic parameters and simulation 


analysis conditions 


参数 数值 
电压 幅 值 /频率 220V/50Hz 
谐 波 3、13 次 谐 波 
白 噪 声 最 大 幅 值 为 1V 
电压 暂 降 (%) 60 
暂 降 起 止 时 刻 /s 0.2 ~ 0.5 
采样 频率 /Hz 3 200 
小 波 函 数 db10 
分 解 尺 度 4 


图 3 为 暂 降 信号 在 小 波 变换 各 尺度 的 传播 特性 。 
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可 以 看 出 ， 暂 降 为 原始 信号 的 一 部 分 ， 分布 在 小 波 
变换 的 各 个 高 频 尺度 上 ， 均 呈现 模 极 大 值 ， 随 着 分 
解 尺 度 的 增加 ， 模 极 大 值 逐渐 增 大 ， 在 高 频段 D, 呈 


现 最 大 值 。 
2 
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图 3 暂 降 信 号 在 小 波 各 尺度 的 传播 特性 
Fig.3 Propagation characteristics of voltage sag at different 
wavelet scales 
图 4 为 噪声 在 小 波 各 尺度 的 传播 特性 。 可 以 看 
出 ， 噪 声 作 为 干扰 信号 随机 分 布 在 小 波 变 换 的 各 个 
高 频 尺 度 上 ， 噪 声 模 值 呈现 负 奇 异 特性 ， 即 在 高 频 
频段 上 较 大 ， 随 尺度 的 增加 逐 新 减 小 ， 在 低频 频段 


上 基本 不 存在 。 
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图 4 噪声 在 小 波 各 尺度 的 传播 特性 
Fig.4 Propagation characteristics of noise at different 


wavelet scales 


图 5 为 谐 波 在 小 波 各 尺度 的 传播 特性 。 不 同 于 
噪声 和 和 暂 降 信 号 ， 各 次 谐 波 分 布 在 指定 高 频频 段 中 。 
3 次 谐 波 主要 分 布 在 Ds 频段 中 ，13 次 谐 波 主要 分 布 
在 DD, 频段 中 ， 基 波 信息 分 布 在 44 频段 中 。 

综 上 所 述 ， 骂 声 、 暂 降 和 谐 波 在 小 波 变 换 中 具 
有 以 下 传播 特性 : 噪声 和 和 暂 降 分 布 于 各 个 高 频 尺度 
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图 5 谐 波 在 小 波 各 尺度 的 传播 特性 
Fig.5 Propagation characteristics of harmonic wave at 
different wavelet scales 

上 ， 而 谐 波 存在 于 小 波 变 换 的 特定 高 次 频段 中 。 如 
图 6 所 示 ， 噪 声 呈 现 与 暂 降 相反 的 奇异 性 ， 对 高 频 
尺度 特别 是 D!/、D;, 等 高 频频 段 具 有 较 大 影响 ， 噪 
声 和 谐 波 (特别 是 高 次 谐 波 ) 的 存在 将 导致 模 极 大 
值 变 弱 ， 影 响 电压 暂 降 定位 。 


暂 降 模 极 大 值 
噪声 模 极 大 值 


Ds 


D, D; 
小 波 频段 
图 6 噪声 与 暂 降 在 小 波 各 尺度 的 模 极 大 值 


Fig.6 The modulus maximum ofnoise and sag at different 


wavelet scales 


3.2 检测 原理 

尽管 噪声 和 谐 波 分 布 在 小 波 变 换 高 频频 段 并 影 
响 极 大 值 的 检测 ， 但 在 最 低频 段 品 声 和 谐 波 儿 乎 没 
影响 。 基 于 此 ， 本 文 提出 了 一 种 利用 双 小 波 变换 
实现 小 波 暂 降 的 实时 检测 方法 ， 在 大 尺度 下 提取 分 
析 电 压 基 波 成 分 ， 在 小 尺度 下 检测 分 析 电 压 突 变 信 
息 ， 从 而 准确 定位 电压 暂 降 的 起 止 时 刻 。 如 图 7 所 
示 ， 双 小 波 检测 的 具体 步骤 如 下 : 

(1) 确定 合适 的 尺度 /， 将 原始 信号 进行 多 尺 
度 小 波 分 解 ， 获 得 各 尺度 的 小 波 系数 。 
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确定 小 波 函 数 WE, 
分 解 尺 度 / 
| 多 尺度 小 波 变换 
i | 获取 各 尺度 小 波 系数 
1 D1 
低频 段 4 单 支 重 构 。 
， | 。 获得 基 波 分 量 
小 流 变 换 2 | 
| 单 尺度 小 波 变换 | 
获得 高 频 分 量 D; 
| 模 极 大 值 检测 
提取 暂 降 特征 信息 


图 7 双 小 波 分 析 检 测 原理 

Fig.7 The principle of double wavelet analysis detection 

(2) 对 最 低频 段 4) 进行 单 支 重 构 ， 得 到 含 暂 降 
信号 的 基 波 分 量 aj。 

(3) 对 基 波 分 量 aj 进行 单 尺 度 小 波 变换 ， 获 得 
高 频 分 量 Di 。 

(4) 通过 检测 高 频 分 量 Di 的 模 极 大 值 ， 确 定 
电压 暂 降 的 起 止 时 刻 。 
3.3 消失 矩 阶 数 

双 小 波 检测 中 的 关键 是 基 波 信息 提取 的 精确 性 
和 和 暂 降 时 刻 检测 的 准确 性 ， 消 失 和 矩 阶 数 表征 小 波 函 
数 有 逼近 信号 的 收敛 率 ， 是 小 波 分 析 的 重要 参数 ， 本 
文选 用 db 系 小 波 函 数 作 为 研究 对 象 ， 分 别 研究 基 波 
信息 提取 的 精确 度 和 和 暂 隆 信息 检测 的 准确 度 与 消失 
和 矩 阶 数 入 的 关系 ， 为 双 小 波 检测 中 小 波 函 数 的 选取 
提供 支撑 。 以 220V 含 30% 电压 暂 降 的 工 频 信 号 为 
例 ， 暂 降 检 测 的 准确 度 和 基 波 检测 精度 分 别 用 偏离 
时 间 A7 和 误差 率 Re 来 表征 ， 即 


AT=|-| (2) 
Rus -< (3) 


式 中 , 大 4 为 检测 的 电压 暂 降 时 间 和 基 波 幅 值 ，、 
4 为 电压 暂 降 发 生 的 时 间 和 基 波 实际 幅 值 。 分 析 结 
果 如 图 8 所 示 。 

从 图 8 可 以 看 到 ， 当 销 失 矩 阶 数 较 小 时 ，db 小 
波 对 基 波 特征 信息 的 提取 误差 较 大 ， 随 着 消失 和 矩 阶 
数 的 增加 ，db 小 波 基 波 特征 信息 的 提取 精度 越 来 越 
高 ， 误 差 越 来 越 小 ， 在 消失 矩 N= 10 时 具有 明显 的 
拐点 ， 另 一 方面 ， 当 消失 和 矩 阶 数 较 小 时 ， 和 暂 降 检测 
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0.025 


0.020 


0.015 


AT/s 


0.010 


0.005 


0.000 


8 消失 和 矩 阶 数 对 误差 率 及 偏离 时 间 的 影响 
Fig.8 The effect of vanishing moments on the error rate and 


time deviation 


的 时 间 误 差 较 小 ， 随 着 消失 和 矩 阶 数 的 增加 ， 暂 降 检 
测 准确 度 逐 渐 降低 ， 二 者 基本 呈现 线性 关系 。 综 合 
考虑 时 间 延 迟 和 提取 精度 影响 ， 本 文采 用 db20 小 波 
进行 基 波 信息 提取 ，db10 小 波 进行 暂 降 信息 检测 。 


4 ”仿真 分 析 


为 了 验证 双 小 波 检测 方法 的 可 行 性 ， 首 先 在 
Matlab 以 第 2 节 中 的 信号 为 例 进 行 分 析 。 根 据 双 小 
波 检测 原理 ， 首 先 以 db20 小 波 进行 多 尺度 分 解 ， 各 
尺度 的 信号 分 量 信息 如 图 9 所 示 ， 在 最 低频 段 4 中 
提取 基 波 信息 分 量 。 然 后 ， 对 基 波 分 量 44 以 db10 
小 波 进行 单 尺 度 分 解 ， 所 得 高 频段 信息 如 图 10 所 
示 ， 通 过 检测 二 次 小 波 变 换 高 频 分 量 的 模 极 大 值 确 
定 电压 暂 降 发 生 的 起 止 时 刻 ， 具 有 较 高 的 准确 度 ， 
能 够 有 效 抑制 噪声 和 谐 波 的 干扰 。 


S 
二 1 1 hn 1 1 1 1 上 1 
.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
A 
-10 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
10 
SS 0 
-10 1 1 1 1 1 1 1 上 1 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
10 T T T T T T Tr T 
SS 0 
-10 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
200 
< 0 
-200 站 


上 1 | 1 1 一 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
l/s 


9 电压 暂 降 信 号 的 四 层 小 波 分 解 
Fig.9 The four-layer DWT of voltage sag signal 
进一步 ， 由 于 小 波 变 换 的 级 联 有 可 能 造成 一 定 
时 间 延 迟 ， 在 Simulink 仿真 系统 下 对 上 述 环境 的 双 
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2001 
-200 : 


0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 


= , , , , ， ， , 1 , 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 
t/s 


10” 双 小 波 分 析 检 测 结 果 


Fig.10 The detect results of double wavelet analysis 


小 波 分 析 方 法 进行 验证 。 双 小 波 检 测 模型 及 环境 如 
11 所 示 ， 暂 降 电源 及 其 干扰 条 件 由 可 编程 序 的 三 
相 电源 提供 ， 双 小 波 分 析 方 法 实现 采用 仿真 系统 自 
带 的 DWT 模块 和 IDWT 模 块 ， 其 中 DWT 模块 实 
现 信号 的 多 分 辩 率 分 解 ，IDWT 模块 实现 单 支 重 构 ， 
选取 A 相 电 压 信 号 进行 双 小 波 分 析 检 测 电压 暂 降 ， 
仿真 结果 如 图 12 所 示 ， 采 用 双 小 波 分 析 方 法 检测 到 


11 双 小 波 变换 仿真 模型 图 


Fig.ll Simulation model of double wavelet transform 


-100 


-200 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 


ts 
(a) 暂 降 信 号 
0.10 
0.05 
SG 0 
-0.05 


-0.10 
0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 0.55 0.60 


ts 
(b) 小 波 变 换 1 高 频 系数 
2 


AQ 0 


0 0.15 0.20 0.25 0.30 0 0.40 0.45 0.50 0.55 0.60 
(c) 小 波 变 换 2 高 频 系 数 
12” 双 小 波 分 析 时 延性 分 析 结 果 
Fig.12 The time-delay analysis results of 


double wavelet analysis 
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的 电压 暂 降 起 止 时 间 分 别 为 0.202s 和 0.502s ,误差 
为 2ms， 具有 和 较 好 的 实时 性 。 


5 结论 


本 文 提出 了 一 种 双 小 波 电压 和 暂 降 检测 方法 ， 首 
先 采用 db20 小 波 对 电压 暂 降 信号 进行 分 解 ， 在 低 
频 尺 度 进行 单 支 重 构 获 取 含 暂 降 的 基 波 信号 ， 再 用 
db10 小 波 对 基 波 信号 进行 奇异 性 检测 ， 精 确 确定 暂 
降 发 生 的 起 止 时 刻 ， 并 通过 Matlab/Simulink 进行 仿 
真 验 证 。 仿 真 结果 表明 ， 双 小 波 检 测 方 法 能 够 有 效 
岂 抑 制 哄 声 和 谐 波 和 干扰， 具有 较 强 的 抗 干扰 能 
在 同时 含有 噪声 和 谐 波 条 件 下 也 可 实现 电压 暂 降 的 
准确 定位 ， 有 效 地 提高 了 小 波 变换 的 检测 性 能 。 
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